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Abstract. Working at twist 2 accuracy and assuming the dominance of the Generalized Parton Distribution H we study the 
helicity-dependent and independent cross sections measured in Hall A, the beam spin asymmetries measured in Hall B at 
Jefferson Laboratory and beam charge, beam spin and target spin asymmetries measured by Hermes. We extract the real and 
imaginary parts of the Compton Form Factor Jf, the latter being obtained with a 20-50 % uncertainty. We pay extra attention 
to the estimation of systematic errors on the extraction of Jff. We discuss our results and compare to other extractions as well 
as to the popular VGG model. 
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The Deeply Virtual Compton Scattering (DVCS) process [1,2] was early recognized as the cleanest way to access 
Generalised Parton Distributions (GPD) and has been so far a very active field of research. Theoretical developments 
are reviewed in Refs. [3, 4, 5, 6] and relevant experimental results are published in Refs. [7, 8, 9, 10, 11, 12, 13, 14, 
15, 16, 17, 18, 19]. 

What we would like is to eventually produce a experimental 3d picture of the nucleon thanks to the knowledge of 
GPDs. Even if getting a precise experimental knowledge of GPDs is a long-term program, we can already perform the 
first fits to get a flavour of the actual sensitivity of observables to GPDs and elaborate robust and efficient extracting 
methods. We address these issues by studying recent JLab data, namely beam spin asymmetries [14] (BSA) and 
helicity-dependent and independent cross sections [12]. These data offer the interesting features of a large kinematic 
coverage and a fine kinematic binning. We also consider Hermes data : beam charge, beam spin and target spin 
asymmetries [15, 20]. 

The first section of this paper describes the evaluation of the ep — > e/?y cross sections and outlines our hypothesis. In 
the second part we explain our fitting strategy. In the third section we discuss the extracted values of Jff and compare 
them to other similar studies. 



Four GPDs H, E, H and E appear at twist 2, but the cross sections depend on the Compton Form Factors (CFF) Jf?, 
Jf? and S. The convention of Ref. [21] is used to define the CFFs. The complex integration kernel yields a real and 
an imaginary part to a CFF. For example the CFF ,¥f satisfies : 



INTRODUCTION 



PRELIMINARY ANALYSIS 



DVCS at leading twist 




(l) 



(2) 



where t, — x B (l + 1 / (2Q 2 )) / (2 — x B + x B t/Q 2 )) is the generalised Bjorken variable [21, 24], x B the standard Bjorken 
variable, Q 2 the virtuality of the initial photon and t the square momentum transfer between initial and final protons. 
3? denotes the principal value of the integral. 



The use of dispersion relations relating real and imaginary parts of a CFF had been discussed by Anikin and Teryaev 
[25, 26, 27], Diehl and Yvanov [28] and Goldstein and Liuti [29, 30, 31]. Dispersions relations are used in the model- 
dependent global fit of Kumericki and Miiller [32] . However the unknown subtraction (the D-term [33]) and the limited 
kinematic range of available data make this constraint rather weak in model-independent fitting. In this work we try to 
minimize the model-dependence of the extraction. Then we consider real and imaginary parts of CFFs as independent 
except when we use an explicit parametrisation of GPDs to compute CFFs. 



ep — > epy observables at twist 2 

The harmonic analysis of ep — > epy cross-sections has so far relied on the 2002 work of Belitsky, Miiller and 
Kirchner [21] (BMK). In this formalism, the interference between the Bethe-Heitler (BH) and DVCS processes was 
treated with a leading order approximation of the BH part. This assumption was removed by Belitsky and Miiller in 
the case of a spinless target [22] and for a (longitudinally polarised or unpolarised) spin 1/2 target [23] and by Guichon 
and Vanderhaegen in the case of a proton target [24]. In all the following, we will use the expressions from the latter. 
At last, the twist 2 approximation is motivated by the claim of early g 2 -scaling observed in Hall A [12] . Moreover the 
^dvcs asymmetry measured by Hermes [20] are small if not zero while they are supposed to vanish exactly at twist 2. 



About current extractions methods 

The methods to extract GPDs from measurements fall today into 4 groups : 

local fits The real and imaginary parts of CFFs are the free parameters. A given kinematic bin (xg,t , Q 2 ) is considered 
independantly of the others. The model dependence is low but some assumptions on the kinematics may be 
necessary to avoid underconstrained problems [38, 39, 40, 41]. It does not allow any extrapolations and in 
particular not those necessary to obtain spatial information in the transverse plane [42] or the quark orbital angular 
momentum in the nucleon [1]. 

global fits CFFs or GPDs have model-dependent expressions, the free parameters of which are fitted. See for example 
[32] for recent fit of DVCS data on unpolarized targets or [43, 44, 45] for recent studies of Deeply Virtual Meson 
Production data. 

hybrid fits This is a combination of the two previous methods, used in [46] and detailed in the present text, 
neural networks fits Work is in progress but preliminary results were recently presented [47]. 



//-dominance hypothesis 

When we consider experiments on unpolarised proton targets we can neglect E, H and E (//-dominance). The 
contribution of these GPDs to helicity-dependent cross-sections is indeed kinematically suppressed [46] and we can 
check (for instance thanks to the VGG model [34, 35, 36, 37]) that Im<§ and /mif are usually smaller than /m^. At 
last for small t and E, we expect H/H to be close to Aq/q i.e. |. 

Thus assuming //-dominance we may hope to extract information on H from BSAs or helicity-dependent cross 
sections with a systematic error of 20 % to 50 this approximation being better at small t. The advantage of this 
approach is the dramatic decrease of the number of degrees of freedom involved in fits. Guidal indeed showed [38] in 
the same kinematic region that it is not possible to extract sensible information about the real and imaginary parts of 
all CFFs by direct fits of helicity-dependent and independent cross sections. 



A direct test of the //-dominance assumption with the VGG model gives an upper bound of 25 %. This is comparable to the typical statistical 
uncertainty on BSAs. 



FITTING STRATEGIES 



The possibility to study GPDs in DVCS rests on factorisation theorems [48] which require a small value of \t\/Q 2 . 
In the following, we restrict ourselves to kinematic configurations for which \t\/Q 2 < 1/2 and justify this choice a 
posteriori. As stated in the first section, we perform two kinds of fits : 

Local fits Re,3>f and ImJf are the free parameters of the fits. 

Global fit We use a "dual model"-like parametrization ; it consists basically in a simultaneous expansion in Gegen- 
bauer and Legendre polynomials. See Ref. [46] for details on the fitting Ansatz. 

We estimate the systematic errors associated to our //-dominance hypothesis by first fitting data setting the subdom- 
inant GPDs to 0, then fitting the same data setting the subdominant GPDs to their VGG value, and computing the 
difference. We estimate the systematic uncertainty related to the global fitting Ansatz by varying the number of poly- 
nomials involved in the expansion. 

A careful fitting strategy is elaborated to handle the high number of free parameters involved in the fits and make 
sure that the numeric implementation is under control. We adopt an iterative (5 steps) fitting procedure. See Ref. [46] 
for full details. 

RESULTS 
Extraction of lm,^ and Re,ff 

The Fig. 1 displays our results for Im,¥f and Re,?? . Both local fits and global fit give results with comparable 
accuracy for Im^P , but as expected the results of the global fits are smoother. This is especially true concerning 
Re.Jtf' : in this case the local fits suffer from large fluctuations of ReJf? with values which fall outside the plot range. 
However, we could not reliably extract values of the CFFs for the larger values of E, with the global fit. 

The results for local and global fits are almost always compatible, which is a strong consistency check. Both rely on 
the assumptions of twist 2 accuracy and of //-dominance. On one hand, local fits suffer from numerical fluctuations 
(the 2-parameter local fits are not constrained enough on some bins) but are almost model-independent. On the other 
hand, global fits are smoother, but suffer from oscillations. That both methods give the same results indicates that 
fluctuations and oscillations are reasonably controlled in the bins for which results are displayed. Since, in both cases, 
the total error bars have the same size, we conclude that our estimation of systematic uncertainties on the global fit 
Ansatz is realistic even if not perfect. 

All fits keep data satisfying \t \/Q 2 < 1 /2. For local fits, changing the maximal value of \t\/Q 2 amounts to dropping 
points. For global fits, the whole results may be changed, but the good agreement between the results of both types of 
fits, and the slow g 2 -evolution of the extracted CFFs, indicate that this restricted kinematic region is suitable for an 
analysis in the GPD framework. 

Our results are explicitly given in Ref. [46]. Error bars are dominated by systematic effects. Typically we obtain a 
relative accuracy of 20 to 50 % on ImJf, which is quite satisfactory under the assumption of //-dominance and given 
the statistical accuracy of JLab data. On the contrary, Re Jf? is still largely undetermined, and is never extracted with a 
precision better than 50 %. 



Discussion 

The Fig. 2 compares our results on Hall A kinematics to twist 2 model-independent extractions [38, 40] of Guidal 
and two extractions with the BMK formalism [12, 32]. Firstly, the use of the GV expressions creates important 
deviations to the extraction of Ref. [12]. Since the extracted combinations of GPDs are not the same, we will not 
make the argument more quantitative. Another discrepancy occur when comparing to the results of the global fit of 
the GPD H of Kumericki and Miiller [32]. Secondly, we obtained results in good agreement with Ref. [38]. Note 
however the surprisingly big shift between the results of Ref. [38] and Ref. [40] induced by the inclusion of target spin 
asymmetries. These points need to be elucidated but we can already mention the following points : 

• Large contributions of H and E are necessary in order to fit JLab Hall A data with the parametrization of 
Kumericki and Miiller. Other DVCS data on unpolarized targets can be included in this global fit under the 




FIGURE 1. g 2 -behaviour ( 1 < Q 2 < 4 GeV 2 ) of the extracted values of Imje (left) and ReJf (right) of local fits and global fit 
on Hall B kinematics : 0.09 < -t < 0.2 GeV 2 (a), 0.2 < -t < 0.4 GeV 2 (b), 0.4 < -t < 0.6 GeV 2 (c), and 0.6 < -t < 1. GeV 2 (d). 
The error bars include both statistics and systematics. ImJt? ranges between and 10 and ReJif ranges between -7.5 and +7.5. The 
black full circles correspond to xb=0. 125, red squares to xg=0. 175, green up triangles to x#=0.250, blue down triangles to x#=0.360 
and magenta open circles to .Y£=0.491. 



//-dominance assumption. 

• The shift between the two extractions of Guidal is attributed to a H contribution. 

• Using VGG and the MIT Bag GPD model [49] we obtained a refined estimate of the systematic uncertainty 
caused by //-dominance. We obtained larger error bars for the extracted CFF on Hall A kinematics, and error 
bars compatible with previous estimates on Hall B kinematics. 

These facts indicate that JLab Hall A data can certainly not be understood under the restricted set of assumptions : twist 
2 + //-dominance. Confirming an explanation involving GPDs other than H will require the inclusion of observables 
more sensitive to these GPDs. Moreover it was already stressed in Ref. [46] that the twist 2 assumption need to be 
checked further. Interestingly the three different methods described in Refs. [32, 38, 40, 46] give results in good 
agreement on Hall B data. The reason of the distinction (discrepancy for Hall A data, agreement for Hall B data) need 
to be elucidated to make progress on GPD extraction in the valence region. Let us remind that JLab Hall A and Hall B 
data are compatible. 



CONCLUSIONS 

Working at leading twist, and assuming //-dominance, we extracted ImJff' and Re,^f with 20-50 % accuracy on Im,3^. 
The good agreement between the results of local and global fits is a strong consistency check. Moreover values of 
Im,^f on Hall B kinematics extracted by three different groups with three different methods are in good agreement. The 
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FIGURE 2. ImJV (up) and ReJP (down) vs t (0.1< -t < 0.4 GeV 2 ) on Hall A kinematics (jt B =0.36 and Q 2 = 2.3 GeV 2 ). ImJf 
ranges between 3 and 8, and ReJif between -5 and 1. We compare our results (left column) to those of Guidal [38] (middle column) 

and Munoz-Camacho et al. [12] (right column). In the latter column M J + jz 1 ^ + j| ~ ^pT^ * s pl otte d, and not J^C. The 
error bars include both statistical and systematic uncertainties. The green (resp. blue) curve is the result of the model-dependent 
global fit of Kumericki and Miiller [32] with (resp. without) The full black square is the result of Ref. [40]. 

comparison of local and global fits and the weak £> 2 -dependency of the results also validates a posteriori the restriction 
to kinematic configurations with \t\/Q 2 < 1/2. Note however that the interpretation of JLab Hall A measurements is 
difficult and probably requires to go beyond the common //-dominance approximation. 

Our results are dominated by systematic uncertainties. The systematic uncertainties related to H-dominance will 
decrease in the future using additional BSA measurements [50] on unpolarised and longitudinally polarised proton 
target which will put stronger constraints on the global fits. These future measurements will also help clarifying the 
interpretation of Hall A helicity-independent cross sections. 
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